Stable propagation of pulsed beams 
in Kerr focusing media with modulated dispersion 
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We propose the modulation of dispersion to prevent collapse of planar pulsed beams which prop- 
agate in Kerr- type self- focusing optical media. As a result, we find a new type of two-dimensional 
spatio-temporal solitons stabilized by dispersion management. We have studied the existence and 
properties of these solitary waves both analytically and numerically. We show that the adequate 
choice of the modulation parameters optimizes the stabilization of the pulse. 



The analysis of the propagation of high power pulsed 
laser beams is among the most active fields of study in 
Nonlinear Optics, where the main dynamics phenomenon 
is the dependence of the refractive index of the materi- 
als with the amplitude of light fields. For propagation in 
materials showing a linear dependence of the refractive 
index with the laser intensity, the mathematical formu- 
lation of the beam dynamics is adequately described by 
the cubic nonlinear Schrodinger equation (NLSE) 0. In 
this case, the excitation of optical solitons is one of the 
most significant phenomena |2|. 

Despite the success of the concept of soliton, these 
structures mostly arise in 1+1-dimensional configura- 
tions. This is mainly due to the well known collapse 
property of the cubic NLSE in multi-dimensional sce- 
narios. This implies that a two-dimensional laser beam 
which propagates in a Kerr-type nonlinear medium, will 
be strongly self-focused to a singularity if the power ex- 
ceeds a threshold critical value, whereas for lower powers 
it will spread as it propagates. Since collapse prevents 
the stability of multidimensional "soliton bullets" in sys- 
tems ruled by the cubic NLSE, a great effort has been de- 
voted to search for systems with stable solitary waves in 
multidimensional configurations 3]. It has been recently 
shown that a modulation of the nonlinearity along the 
propagation direction in the optical material can be used 
to prevent the collapse of two-dimensional laser beams 
Q . The concept has been extended to the case of several 
incoherent optical beams p and to the case of matter 
waves 0. 

In this paper we use a similar idea to stabilize against 



collapse pulsed laser beams propagating in planar waveg- 
uides. This configuration is of great importance, as it 
is the common experimental procedure used for exciting 
spatial optical solitons. Instead of making a modulation 
of the nonlinearity our idea is to act on the chromatic dis- 
persion term of the NLSE. Thus, our procedure resembles 
that used to obtain dispersion- managed temporal solitons 
in optical fibers. 

We consider the paraxial propagation along z of a 
pulsed beam of finite size in time t and spatially con- 
fined by a waveguide along the ?/-axis. Thus, diffraction 
only acts in the x direction and is balanced by the self- 
focusing nonlinearity given by a refractive index of the 
form: n = no + ri2\E\ 2 . The dynamics of the slowly- 
varying amplitude of the pulse E is described by a 1+2D 
NLSE of the form: 

ko ox z no 

a) 

where ko,k and k f are respectively the wavenumber in 
vacuum, the inverse of the group velocity, and the group 
velocity dispersion coefficient. 

As we mentioned above, for continuous beams an ad- 
equate modulation of n2 along z prevents collapse of the 
light distribution and yields to stabilized two-dimensional 
solitons. From a mathematical point ow view an equiv- 
alent effect can be achieved by suitably modifying the 
second order derivatives of the propagation equation • 
However, it is physically impossible to change the sign 
of diffraction of a continuous beam during its evolution. 
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Nevertheless, for the case of pulses, the temporal deriva- 
tive in Eq. JIJ can be modulated by changing the dis- 
persion parameter fcg, which can take positive as well 
as negative values. In this paper we will show that it 
is possible to stabilize a pulsed beam against collapse by 
modulating only the dispersion during beam propagation, 
without altering diffraction. 

For a linearly polarized Gaussian pulsed beam with 
power P = J \E\ 2 dzdx, an initial beam waist wo yX and 
temporal width wo,t, it is useful to write the above Eq. 
(PJ in adimensional form by means of the Fresnel length 
F = koWQ x in the "reduced time" frame. Thus, we 
make the rescaling: r = (t — k' z)/wo yt , n = z/F, x — 
x/(F/k )^ 2 and u = (w^ x w^ t /k^) 1/2 E / 'P 1 / 2 . Then, if 
we consider that dispersion is modulated by a periodic 
function along propagation, we will finally yield to a gen- 
eralized NLSE of the following form: 



.du 1 .d 2 u 
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+ g\u\ 2 u = 0. 
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Here, d{rj) = Fk^/w^ and g = kon2PFk f /n wo jX wo jt . 
It is not obvious what kind of periodic function g(rj) will 
be more adequate to stabilize pulsed beams. A simple 
choice (and very natural for applications) is to take d(rj) 
as a piecewise constant function of the form d{rj) = d a 
for n < 2r] a and d(rj) = dt for 2n a < n < 2775 (see top of 
Fig. ^) where d a , db, r] a and 775 are free parameters that 
can be varied in order to fit experimental requirements 
and to optimize the stabilization of the pulsed beam. 

To get some insight on the dynamics given by Eq. 
we have performed an analytic study by means of the 
time-dependent variational approach |s[ . Notice that Eq. 
(J2J can be obtained from the Lagrangian density: 
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2C = i (uu* — u*u) + d{rj) 



where dots denote derivative with respect to rj. 
choose a Gaussian ansatz of the form: 
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The n— dependent parameters in the above equation have 
the following meaning: A is the amplitude, w x ,w T are 
the spatial and temporal widths and /3 X , (3 T are the initial 
curvature and chirp. Although Gaussians are not exact 
solutions of Eq. our choice simplifies the calculations 
and is a usual lithg distribution in experiments. The 
standard variational calculations |8| lead to the equations 
describing the dynamics of the pulsed beam spatial and 
temporal widths: 



VJy = 
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w x 
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(5a) 
(5b) 



Numerical simulations of Eqs. ((5J) show that the beam 
can be stabilized for many choices of the model param- 
eters. These equations also predict two types of oscil- 
lations for width w T : a fast one due to the modulation 
of dispersion and a slow one which almost coincides with 
the variation of w x . This low frequency oscillation is gen- 
erated by the internal nonlinear dynamics of the system. 
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FIG. 1: Top: Sketch of the propagation of a pulsed beam 
through a system with periodic modulation of dispersion 
made by joining two different planar Kerr waveguides. The 
forward- backward trajectory can be obtained by reflection at 
mirrors Ml and M2. Bottom: oscillations of the peak am- 
plitude of stabilized pulsed beams propagating in the above 
system. In (a) the upper curve corresponds to a value 
g = 1.074g c (surface plots in top figure). For the lower curve 
g — 1.034g c . The modulation is of the form: d a = 1 + 8, 
d b = 1 — 8, 2r\ a = 2r\ h = 0.04. In (b) and (c) details of the 
fast oscillations are shown. 



Variational models usually provide only simple and in- 
tuitive predictions that must be complemented with the 
integration of Eq. (0). All the results to be presented in 
this paper are based on direct simulations of Eq. us- 
ing a split-step Fourier method on a 520 x 520 grid and 
absorbing boundary conditions to get rid of radiation. 
In Fig. H we plot the evolution of the peak amplitude of 
two different pulsed beams in a system sketched in Fig. 
ID (top). In Fig. [IJa) the upper curve corresponds to a 
value g = 1.074g c , being g c > the critical value of the 
nonlinearity for collapse. For the lower curve g = 1.034g c . 
The modulation is of the form: d a = 1 + 8, d^ = 1 — 8, 
2r] a = 2^5 = 0.04. In (b) and (c) details of the fast 
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oscillations are shown. 

In Fig. [3 we plot the width oscillations of the sta- 
bilized soliton from Fig. ^ with g — 1.074g c . As the 
variational equations predict, the spatial width only dis- 
plays a low frequency oscillation, whereas the temporal 
width presents small amplitude fast variations over the 
same low frequency main oscillation. This is clearly seen 
in the frequency spectra displayed in Fig. Efd) and Fig. 
He). 
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FIG. 2: (a) Oscillations of the spatial (lower) and temporal 
(upper) widths of the stabilized soliton with g = 1.074g c from 
Fig. 0^a). Insets (b) and (c) display details of the curves. 
Insets (d) and (e) show the frequency spectra of the spatial 
and temporal width oscillations, respectively. 

In Fig. |3| we plot the results of simulation of Eq. © 
for the same input pulsed beams of Fig. ^ Now we 
take waveguides of half size, resulting in a frequency of 
modulation which is twice that of Fig. ^ Stabilization 
of both beams is achieved, showing the robustness of this 
idea. 

As a conclusion we can say that in this paper we 
have described two-dimensional spatio-temporal solitons 
which are stabilized against collapse by means of disper- 
sion management. We have studied their stability and 
properties both analytically and numerically. We have 
shown that the adequate choice of the modulation pa- 
rameters will optimize the stabilization of the pulse. Our 
results are of importance in the field of high power pulse 
propagation in nonlinear optical materials. 

A remaining open question concerns the stability of 
higher-dimensional beams. In fact, at the present time it 
is not clear the mechanism of stabilization of fully three 
dimensional solitons which have been predicted to col- 
lapse in the case of modulation of nonlinearity p. Thus, 
our results on dispersion management would be an im- 
portant step forward in stabilization of light bullets. 




FIG. 3: Oscillations of the peak amplitudes (a) and widths 
(b) of the same input solitons from Fig. In this case the 
modulation is given by d a = 1 + 8, db = 1 — 8, 2r] a — 2r]b — 
0.02. In (b) the oscillations of the spatial and temporal widths 
are overlapped (insets). 
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